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A double peak structure in the magneto-phonon resonance (MPR) spectrum of uniaxial strained 
graphene, under crossed electric and magnetic fields, is predicted. We focus on the F point optical 
phonon modes coupled to the inter-Landau level transitions 0 ±1 where MPR is expected to be 

more pronounced at high magnetic held. We derive the frequency shifts and the broadenings of 
the longitudinal (LO) and transverse (TO) optical phonon modes taking into account the effect of 
the strain modihed electronic spectrum on the electron-phonon coupling. We show that the MPR 
line for a given phonon mode acquires a double peak structure originating from the two-fold valley 
degeneracy lifting. The latter is due to the different Landau level spacings in the two Dirac valleys 
resulting from the simultaneous action of the inplane electric held and the strain induced Dirac cone 
tilt. We discuss the role of some key parameters such as disorder, strain, doping and electric held 
amplitude on the emergence of the double peak structure. 
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I. INTRODUCTION 


The discovery of the anomalous quantum Hall effect 
(QHE) in graphene^ was the fingerprint of the presence of 
two-dimensional (2D) relativistic-like electrons, so called 
Dirac electrons, which are at the origin of the striking 
properties of graphene. 

In the presence of a magnetic field perpendicular to the 
graphene sheet, the electron spectrum transforms into 
discrete Landau levels (LL) with typical energy e„ = 


sgn{n) y/2\n\ where n is an integer, vp is the Fermi 


velocity, Ib = 



is the magnetic length and sgn(n) = 


-!-(—) in the conduction (valence) band^’^. 

To study the possible inter-Landau level excitations, 
magneto-Raman spectroscopy has proven to be a good 
probe which is remarkably sensitive to the electron- 
phonon coupling^’^. Such coupling in graphene has been, 
recently, a hot topic of study^. The case of the center- 
zone doubly degenerate E 2 g optical phonon modes and 
the corresponding Raman G band are found to be of 
particular interest^’®^®. By tuning the magnetic field, 
the G band shows a rich structure characterized by an¬ 
ticrossings known as magneto-phonon resonance (MPR) 
effect®The latter occurs when the energy between 
two LL matches that of the G band phonon. The anti¬ 
crossing behavior of MPR was first predicted by Ando® 
who found that the frequency shifts and the broadenings 
of the P point optical phonon modes, in graphene, are 
marked by a singular behavior which can be tuned by 
doping. Goerbig et al}^ proposed that the fine structure 
resulting from the anticrossings of the MPR is expected 
to split due to distinguishable behavior of circularly po¬ 
larized phonons. The properties of the split fine structure 
are found to be hlling factor dependent^®. This depen¬ 
dence was proposed as a probe to measure the strength of 


the electron-phonon coupling and to resolve the phonon 
polarization^®. 

Several magneto-Raman measurements have given ev¬ 
idence of the predicted fine structure properties of the 
MPR of the graphene G band®dH7 particu¬ 

lar its filling factor and light polarization dependences. 
The most pronounced effect of MPR was found, at high 
magnetic field, for the inter-LL transitions 0—^1 and 
— 1 —>• 0^®’^®’^^. A fine structure of the G band was also 
observed in non-resonant regime^®. 

Recently, the filling factor dependence of the MPR was 
studied by controlling the position of the Fermi level in 
gated graphene^®’®®. 

Some puzzling features of the measured MPR struc¬ 
tures remain unsolved^®’^^. In particular, the appear¬ 
ance of an additional component in the MPR structure^^ 
is still not clearly understood. Spatially inhomogeneous 
carrier densities and strain were proposed to be at the 
origin of this component 

The behavior of the MPR in graphene under uniax¬ 
ial strain was theoretically studied in Ref. 21 taking into 
account the effect of the strain on the optical phonon 
spectrum. The interplay between strain and MPR gives 
rise to a Raman line splitting and a light polarization 
preferences®^. In Ref. 21, the authors did not consider 
the strain-induced modification of the electron-phonon 
coupling. They argued that, considering a linear elec¬ 
tronic spectrum with strain induced anisotropic Dirac 
velocities, may simply result in a weak renormalization 
of the magnetic field®^. 

However, it has been shown that the electron-phonon 
interaction is sensitive to the modified electronic spec¬ 
trum due to the strain®®. In particular, the observed 
light polarization dependence of the G band in strained 
graphene was ascribed to the strain induced anisotropy 
of the electronic dispersion®®’®®. 
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The natural question which arises is how the MPR of 
the center zone optical phonon modes is modified if one 
takes into account the effect of the strain on the electron- 
phonon interaction. 

Under strain, the Dirac valleys of the honeycomb lat¬ 
tice are no longer located at the K and K' points of the 
Brillouin zone^. They are shifted away, depending on the 
strain amplitude. Moreover, the Dirac cones, which are 
isotropic in undeformed graphene, become tilted. 

A striking valley dependence of LL energy in graphene 
under crossed electric and magnetic fields was predicted 
by Lukose et The authors showed that, for a uni¬ 

form electric field applied along the graphene sheet, the 
spacing between LL decreases and may eventually vanish 
leading to a collapse of LL by increasing the electric field 
amplitude. This behavior is substantially different from 
that expected for a non-relativistic two dimensional elec¬ 
tron gas where the LL spacing is found to be independent 
of the electric filed^^. 

Goerbig et found that crossed electric and mag¬ 
netic fields can induce a lifting of the twofold valley 
degeneracy of LL in tilted Dirac cones systems as in 
the organic conductor a-(BEDT-TTF) 2 l 3 where BEDT- 
TTF stands for the organic molecule bis(ethylenedithio)- 
tetrathiafulvalene. This salt is a layered material ex¬ 
hibiting Dirac-like electrons under pressure or uniaxial 
strain^®. 

Goerbig et al?^ showed that, in the presence of a 
transverse magnetic field, the inplane electric field gives 
rise to an additional effective tilt term which is valley 
dependent. The authors proposed to probe such valley 
filtering effect by infrared transmission spectroscopy. 

Valley degeneracy lifting was recently investigated in 
graphene. Several ways were proposed to control the 
valley polarization^^ such as introduction of topological 
line defects^® and strain engineering^®. The latter 
is expected to generate pseudomagnetic fields with 
opposite sign in the two valleys. It is worth noting 
that the theory of the strain induced pseudomagnetic 
field and the underlying approximations were recently 
discussed in Ref. 30. 

A strained graphene barrier with artificially induced 
electron mass and spin-orbit coupling (SOG) was pro¬ 
posed as a possible spin-valley filter®^. The electron mass 
and the SOG are expected to result from adatoms or sub¬ 
strate deposition®^. 

The valleytronics in Dirac electron-like systems has 
emerged as a possible way to encode carrier information 
as in spintronics. One needs to remove the valley degen¬ 
eracy to filter the electron of one valley^^. 

In bismuth, in which the Dirac Hamiltonian was 
first introduced to describe the “Dirac electrons” in 
a condensed matter®^, valley degeneracy among three 
equivalent electron pockets was spontaneously removed 
at low temperature and/or high magnetic field. Loss of 
threefold symmetry was reported both in the magneto¬ 


transport®® and in the thermodynamic properties®^. 

In this paper, we consider a honeycomb lattice under 
uniaxial strain in the presence of crossed inplane elec¬ 
tric field and a transverse magnetic field. We investigate 
the behavior of the MPR spectrum taking into account 
the strain induced dependence of the electron-phonon in¬ 
teraction. We focus on the center zone longitudinal op¬ 
tical (LO) and inplane transverse optical (TO) phonon 
modes corresponding to the G band in graphene. We 
consider the coupling between the optical phonon modes 
and the electron-hole pairs resulting from inter-LL tran¬ 
sitions 0 ^ ±1, where the effect of the MPR is more 
pronounced®®. The main result of this work is the pos¬ 
sible occurrence of an electric field induced double peak 
structure in the MPR spectrum. This structure is the 
signature of the valley degeneracy lifting. We also found 
that, considering strain modified electron spectrum, does 
not simply yield, as argued in Ref. 21, to a weak renor¬ 
malization of the magnetic field but contributes to the 
fine structure of the phonon lines induced by the strain. 

The paper is organized as follows. In section II we 
derive the phonon self-energy in tilted Dirac cone systems 
under crossed electric and magnetic fields. The shifts 
and the broadenings of the MPR modes are discussed in 
section III. A brief summary is given in section IV. 


II. OPTICAL PHONON SELF-ENERGY 

It has been argued that the honeycomb lattice under 
uniaxial strain could be described by the minimal 
form of the generalized 2D Weyl Hamiltonian®’®®^®® 
characterized by anisotropic electronic velocities and a 
parameter measuring the tilt of the Dirac cones. The tilt 
effect is not significant in deformed graphene contrary 
to the case of the organic conductor q;-(BEDT-TTF) 2 I 3 . 

The 2D Weyl Hamiltonian can be written as®’®®’®®": 

H^{k) = ^ (^WQ.ka^ + WxkxCT^^ + Wykya^ ( 1 ) 

where ^ = ± is the valley index, k = {kx,ky) is the 
momentum vector, wq = {wox,woy)j cr® = 1, and 
are the 2x2 Pauli matrices, wq is responsible of the 
tilt of Dirac cones whereas Wx and Wy give rise to the 
anisotropy of the energy dispersion. 

Gonsidering a uniaxial strain along the y direction, the 
distance between first neighboring atoms along the y axis 
changes from a to a' = a -|- <5a = a(l -|- e), where e = ^ 
is the lattice deformation which measures the strain am¬ 
plitude. e is negative (positive) for compressive (tensile) 
deformation. 

The lattice basis changes from (di = y/3aex,a2 = 
^acx + \aey) to (di = 'Jiaex,a 2 = ^acx + 
a (I -I- e) By) which corresponds to a quinoid lattice®®. 
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In this case, the parameters of the generalized Weyl 
Hamiltonian (Eq. 1) take the form 

Wx = \/iat sin 9, Wy = —ta cos 9 1'a{l + e) 

Wox = 2V5a{tnnn sin 20 + sin 0 ), and woy = 0 ( 2 ) 

where 9 = arccos^— t and tnnn denote the hop¬ 
ping integrals, respectively, to the first and to the sec¬ 
ond neighboring atoms along the bond directions which 
are not affected by the strain, t' is the hop¬ 

ping integral between first (second) neighbor atoms sep¬ 
arated by the deformed lattice distance a' = a(l -I- e) 
(||a 2 || = V3a{l + |e)). t’ and t(j„„can be written as^ 


It is worth noting that we only consider the effect of a 
uniaxial strain and not a shear one, to avoid complexity 
related to the contribution of Poisson ratio. The 
latter could be neglected as a first approximation^^’^®. 
Moreover, we did not take into account the effect of 
the strain on the phonon dispersion to clearly identify 
the role of the electron-phonon interaction in the strain 
induced MPR fine structure. These approximations were 
discussed in Ref. 22. 


A. Landau levels in tilted Dirac cones at zero 
electric field 


t' =t+ ^Sa 
da 


t' = t 

^nnn ^ri 


dtr, 


da 


-5a 


( 3 ) 


Given the Harrison’s law^ ^ , t' becomes : t' = 

oa a ’ 


t(l- 2 e). 

The electronic spectrum, corresponding to the Hamil¬ 
tonian, given by Eq. 1 is : 


£\{k) = wo-k -f X^Jwlkl + wlkl (4) 


In the presence of a magnetic field transverse to the 
lattice plane B = Be^ and in the Landau gauge A = 
(0, Bx, 0), the electronic spectrum splits into Landau lev¬ 
els (LL) given by:"^*^ 

,/WrW,, , - 

En = sgn{n)^ -\/273|n| (6) 

,— Lb _ 

where = \/^ is the magnetic length, 7 = ^/l — 

and tco is the tilt parameter given by: 

]j \W^ J \Wy J 


where A = ±1 is the band index. 

For small strain amplitude, Wx, Wy and wqx could be 
approximated by^®: 


The corresponding eigenfunctions are of the form^^o 

{f]nX) = Fnx{f) = -^= exp{i^)(j)i{x) ( 8 ) 
V^y '■B 


w. 


3 

^-at 






3 

^-at 


^ 4 


Wqx - 0-Gwxe 


where X = kyl^ is the center of coordinates and Ly is 
the crystal length along the strain direction y. are 

(5) given by'*°: 

I 


= 


2 ^ 1 +^ 


^ + 7 


f\n\ 


sgn(n) f {1 + 57 )e"*^'^* 


2 ^ 1 +^ 


-Wo 


f\n\ — l 


4 = 


^/W+W) V ^ + 7 


/o 


(9) 

( 10 ) 


where tan tit = . In the present case, dt = 0 since 

' WyWQx 

woy = 0 (Eq. 2). 

The f\n\ functions are given by^° 



It is worth stressing, that we have derived Eqs. 9 and 10 
using the method proposed by Himura et However, 
the form of function (Eq. 9) is slightly different from 
that found in Ref. 40 where a sgn(n) was missing and a 
factor of ^ -I- 7 should be replaced by 1 -I- ^ 7 . 


H\n\ being the Hermite polynomials and u 
^{x- l%ky) - ^sgn{n)wolB\/^\ [Refs. 24, 40]. 


( 11 ) 
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B. Landau levels in tilted Dirac cones under a 
uniform electric field 

The electron spectrum in tilted Dirac cone systems, 
in the presence of an inplane electric field E = Ecy is 
characterized by an electric field dependent tilt vector^® 

E x. B 

W^iE) = W^y) =Wo- g2~ (12) 


The spacing between LL is then valley dependent: it is 
larger in ^ = + (C = —) valley for a tensile (compressive) 
deformation as we shall show in the next section^®. 

The LL eigenfunctions (t>n{^) have the same form as 
those given by Eqs. 9 and 10 but with valley dependent 
tilt parameters:^® = y^l — w|(E) where w^i^E) is 

given by Eq. 13. 

C. Electron-phonon interaction Hamiltonian 


with wq = {wqx, woy = 0) is the tilt vector in the absence 
of the electric field (Eq. 2). 

The tilt parameter is given by^®: 


w^{E) = 



(13) 


The tilt angle under electric field is written as^® 
tan$«(E) = 

Under a uniform electric field E = Ecy, the tilt pa¬ 
rameter, given by Eq. 7, becomes: 


Wi{E) = (14) 

Wx ^ 

The LL spectrum is given by^"'’’^® 

4,{E) = sgn{n)^EEl [1 _ w|(E)] ^ 

(15) 


The electron-phonon Hamiltonian in graphene can be 
written as"*^^: 



(16) 


where M is the mass of the carbon atom, (^f./x) is 
the creation (annihilation) operator of phonon with wave 
vector q = (qy^^qy). /3 = = -f|i, is the 

r point optical phonon frequency of the mode /r in the 
absence of electron-phonon interaction, fj, = L = T) 
for LO (TO) phonon mode. 

In undeformed graphene ujqt = ujql = This 

degeneracy is lifted in the strained graphene due to the 
symmetry breaking. However, we assume, for simplicity, 
that ujqt — woL — wq. This turns out to neglect the 
effect of the strain on the phonon dispersion since we 
are interested in the role of the electronic dispersion on 
the magnetophonons^^. 

In Eq. 16, the matrices V^{q) near the Dirac points D 
(^ = -b) and D' (^ = —) are given by^^.: 


Vr^i^ 


WxW, 


y 1 +^ q ;' cos (/3(^ 




0 


^ cos^(g) _|_ 

0 


(17) 


where is the phonon angle tan ip(q) = w'y = 
Wy — 2et’a{l -b e) ~ Wy{l — |e) and a' = Jwx/w'y. 


D. Optical phonon self-energy 

The total self-energy H^ of a given optical phonon 
mode (/r =LO, TO) is due to the interaction of phonons 


with electrons in both valleys (^ = ±), it is written as® 

= E ( 18 ) 

The r point optical phonon self-energy of a mode /r, 
due to interactions with electron-hole pairs near the val¬ 
ley is given by® 
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0,w) = -gs 


hS' 

NMwq 



n,n' X 


f {ei{E)) - f (ei,{E)) 
tvjj + £^, {E) - ei{E) + iS 


(19) 


where we considered the valley dependent LL spectrum 
under electric field E = Ecy (Eq. 15). The functions 
(r|nX) are the LL eigenfunctions given by Eq. 8, f{x) 
is the Eermi-Dirac distribution, gs is the spin degeneracy 
and S = ^, T being the scattering time. 

The sum over X gives rise to the degeneracy of the 

o/ 

LL Ol = T^^TT, where S is the surface of the strained 
graphene which can be expressed in term of the unde¬ 
formed surface S as S' = iV||ai x 02 !! ~ ^ (l -I- |e). 

The self-energy Eq. 19 exhibits a resonance behavior if 
the following condition is fulfilled 

hu;o + ei,{E)-ei{E):^0 ( 20 ) 

where Hljq = 0.196 eV is the bare optical phonon energy. 

Since the E phonon modes induce a vertical inter-LL 
transitions, the momentum dependent term in the 

energy spectrum Eq. 15, responsible of the inclination of 


the LL, can be omitted. The inclination effect is not rel¬ 
evant since, in a given valley, the LL spacing is constant. 
Moreover, the contribution of such effect compared to the 
resonance energy can be estimated by calculating the ra- 
tio where Ae = hkQ A, /cq being the momentum value 
corresponding to a resonant inter-LL 0 ^ ±1 transitions: 
fiwo hvpko where vp ~ 10® m.s“^ is the average Fermi 
velocity in graphene. We then obtain 


hujQ vpB 

The MPR associated to inter-LL 0 ^ ±1 transitions, 
take place around B= 30 Pqj. electric field 

of the order of E ~ 10® V/m, we have ~ 3%. The 
inclination of the LL can then be neglected in this case. 

Using Eqs. 9 and 10, the matrix elements in Eq. 19 
take the form: 


\{nX\V[{cj)\n X)f = w^^WyK ^x 

I (“^) + + {a' cosg}f [(^ + qH-E^))^ - w|(£;)] ^ | 

\{nX\V^{^\n'X)\'^ = w^w'yK^K^^ [V'l.UI-i + x 


here = ip{q} is the phonon angle, ^ As mentioned in Refs. 9, 41, the contribution of w = 

V l .1 Q should be subtracted from the phonon self-energy to 


and = 




avoid a double counting problem. We then obtain: 


Uliq^6,uj) = -CAl E 

A,A'=±1 n=0 


f {XeijE)) - f (A'£i+i(i:;))] (a4(E) - X'ej^.iE)) 


1-AA' 


{hu; + iS)^ — (^Xen{E) — X'e^^^{E)^ ^n+ii^) + £n{E) 


(23) 


QTiri — 1 mUo constant C is given by: 

'-'n^O ~ (^4(1-|-{7«(£;)) ) '-'n=0 — 

C = 5.36^3 (^0 

~ gs^^Vs (^ 


h 1 

2Ma^uJo huiQ 

h 1 

2 J 2Ma^uJo fiLOo ’ 


( 24 ) 
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here oq = a-\/3 is the lattice parameter in the undeformed 

I 


system. The coefficients are given by: 


A{ = wWy 


Sin if 


[(C + 7^(^^))^ + + (a'cos^)^ [(C + 7^(£^))^ - 


A^j, = w^w'y I [{^ + + (a'sint/j)^ [{( + - wl{E)y 


(25) 


Since these coefficients are valley dependent, the 
phonon frequency shifts and the corresponding 

broadenings are also expected to depend on the val¬ 
leys and they are given by: 

Aw^ = - Wo = (g, Wo) 

r« =-i3n«(g,wo) (26) 

Let us, first, focus on the zero electric field case. The 
resonance occurs, according to Eq. 23, when the phonon 
frequency is equal to the energy separation between two 
LL: Hloq = Ae„ — A'e„_|_i. 

For interband (A = — A') and intraband (A = A') transi¬ 
tions between LL, the resonance takes place, respectively, 
at magnetic energy huf^ and /iw^” satisfying®: 

= y/WxWy^^^^^— [v^nTT ± y/n] 

Lb 

= "" ^7 ^ [Vn-I-1 ± Vn] (27) 

6 at 

where fvjjB = [Ref. 9]. 

In the case of uniaxial strained graphene, hjj^ could 
be written, in the limit of small strain amplitude, as: 

Sw^* ~ huB ^1- ^1- ^(0.6e)^^) [\/n -h 1 ± y/n] 

(28) 

where we considered the approximate values of Wx-, Wy 
and wq at zero electric field (Eq. 5). 

Eq. 28 shows that the resonance takes place at higher 
(lower) magnetic field for tensile (compressive) defor¬ 
mation compared to the undeformed case. According 
to Eq. 28, and disregarding the strain modified phonon 
spectrum, the electron-phonon interaction in strained 
graphene induces a shift in the anticrossing position of 
the magnetophonon spectrum. This shift is due to the 
strain modihed electron-spectrum. The shifted position 
of the MPR, compared to that in the undeformed lattice, 
is an indication of the presence of deformation. The 
latter may accidentally occur during the fabrication pro¬ 
cedure of the graphene sheet by exfoliation or chemical 
vapor deposition. 


Let us now turn to the possible evidence of the valley 
filtering in the MPR spectrum. According to Eqs. 23 and 
25, the frequency shifts of a given phonon mode is valley 
dependent. However, at zero electric field, the resonance 
condition, which is responsible of the anticrossing struc¬ 
ture of the MPR, does not depend on the Dirac valleys 
(Eq. 27). One should, then, not expect in the absence of 
the electric held, any signature of the valley degeneracy 
lifting in the MPR spectrum. 

However, under an inplane electric held, the resonance 
condition is valley dependent (Eq. 20). The self-energy is 
then expected to have a singular behavior at two critical 
magnetic energy and corresponding to the 

valleys ^ = ±. If the LL spacing in the ^ = -I- valley is 
larger than in the ^ = — (Eq. 15), the resonance due to 
inter-LL transitions in ^ = -|- valley will then appear be¬ 
fore that associated to ^ = —, by increasing the magnetic 
held, since huj^ is given by: 

huio = [l - w|(A)] [y/n + l + 

O CLl 

(29) 

A double peak in the MPR lines may, then, appear in 
the presence of the inplane electric held. 

It is worth stressing that, besides the phonon line split¬ 
ting induced by valley degeneracy lifting, one may ex¬ 
pect a change in the MPR spectrum due to the elec¬ 
tron tunneling between valleys. This effect takes place 
at high magnetic held or as the merging of Dirac cones 
is approached^^. However, such effect is not relevant 
in graphene since a large strain amplitude is required 
(~ 30%) to observe it. 

Due to the electron intervalley tunneling the zero LL 
splits into two levels 0 ®' and 0 “ separated by an energy 

gap42.43 

AE oc exp 4-\/2-g-^ (30) 

where a=. — cB.l + e and H is a dimensionless pa- 

y '^y 

rameter given hy B = ^ with Hg = ^ ■ 

The intra-LL transitions O®' —> 1 and —1 —>■ 0“ 
will appear at different magnetic helds compared to the 
inter-LL transitions 0“ —> I and —1 —O®'. In the 

















absence of electron intervalley tunneling, the —1 —> 0 
transition is forbidden for a totally filled n = 0 level. 
The corresponding Raman line obtained with circularly 
polarized light disappears and only the line ascribed to 
0 —1 transition persists at high magnetic field. 

In the presence of the intervalley tunneling, besides 
the disappearance of the —1 —> 0 line, the intensity of 
the inter-LL 0“ —> 1 transition will be reduced since 
the intra-LL transition 0+ —> 1 will no more contribute 
to the inter-LL line 0“ —>■ 1. The O’*' —>■ 1 intra-LL 
line may emerge at higher field compared to the inter-LL. 

In graphene, the distance between first neighboring 
atoms is a = 1.42 A. For a critical compressive strain 
e = 30%, before sample cracking, the anisotropy param¬ 
eter is a ~ 0.7, which gives rise to Bq ~ 10"^ T (Eq. 30). 
A magnetic field of 10^ T is, then, needed to detect the 
change in the MPR spectrum of graphene due the val¬ 
ley electron tunneling. However, in the organic salt a- 
(BEDT-TTF) 2 l 3 , the anisotropy parameter a could be of 
the order of 0.4 and Bq ~ 660 T (a ~ 10 A). The changes 
in the MPR spectrum induced by the electron valley tun¬ 
neling may then be observed in a-(BEDT-TTF) 2 l 3 at a 
field amplitude of H > 15 T. 


III. RESULTS AND DISCUSSION 

We first focus on the case of uniaxial strained graphene 
in the absence of the electric field and then turn to the 
study of the electric field induced valley degeneracy lift¬ 
ing. 


A. MPR in uniaxial strained graphene at zero 
electric field 

We numerically derive the frequency shifts and the 
broadenings of the LO and the TO modes correspond¬ 
ing to 0 ^ ±1 transition (Eqs. 23 and 26). 

We consider linearly polarized phonon modes along 
(LO) and perpendicular (TO) to the phonon momentum 
direction {cos ip{q), sin ip{q)) where ip{q) is the phonon 
angle. We denote hereafter, (p{q} = </). For ^ the 

LO mode is along the strain direction. 

Figure 1 shows the frequency shifts of the LO and the 
TO modes in undeformed lattice and for a tensile defor¬ 
mation of 5%. The calculations are done in the undoped 
regime. 

As shown in Fig.l, the strain splits the degenerate 
phonon line of the undeformed lattice into two lines cor¬ 
responding to the two linearly polarized LO and TO 
modes. Kashuba and Fal’ko^^ showed that the Raman 
signal splits under strain but did not address the strain 
dependence of the anticrossing position. 

According to Fig.l, the anticrossings marking the 
MPR, corresponding to the 0 ^ ±1, take place at 



FIG. 1. Frequency shifts of LO and TO modes as a function 
of the magnetic energy Aus in undoped graphene {Ep = 0) 
for a disorder amount of = 0.02 and under a tensile strain 

TUJQ 

e = ^ = 5%. The LO mode is along the strain axis. The 
thin solid line is the result for the undeformed case. 

slightly shifted magnetic field value compared to the 
undeformed case where MPR occurs at Hujb = huiQ- The 
shift is due to the strain modified cyclotron energy. 

Indeed, under uniaxial strain, the resonance condition 
of the 0 ±I transition in graphene is of the form 
(Eqs. 27 and 28) 

-n-l- 2 ^/WxWy , 3 

^hwB (l-^(0.6e)2^(3I) 

The magnetic field at which the anticrossing occurs is 
then huJB ~ ^0 (l + (l + |(0-6e)^) which is en¬ 
hanced (reduced) for tensile (compressive) deformation. 

In undeformed graphene, the resonance corresponds to 
Bq = 30.1 T®. A strain of 5% is expected to shift this 
value to Bq ^ i?o (l + gc) ~ 30.6 T. It turns out that a 
shift of the resonant magnetic field, at which the MPR 
occurs, could be a signature of the presence of strain in 
the sample. 

Numerical calculations show that the splitting of the 
phonon line induced by strain is not significant for small 
strain values. 

For clarity, we consider in the following strain values 
|e| > 10% as shown in figure 2 where we plot the fre¬ 
quency shifts of the TO mode in the undoped regime for 
different strain values. The higher the strain, the larger 
the shift of the anticrossing position. 

The anti-crossing structure characterizing the phonon 
frequency line (Figs. 1 and 2) is due to the resonance 
condition giving rise to the singular behavior of the 









FIG. 2. Frequency shifts of the TO mode as a function of 
the magnetic energy huJB for different strain values in the 
undoped graphene and for a disorder amount 5 = 0.02 fejo- 
The TO mode is transverse to the strain axis. The thin solid 
line is the result for the undeformed case. 


phonon self-energy (Eq. 19). According to this equation, 
the resonance condition depends on the disorder param¬ 
eter S. One should then expect a disorder dependence of 
the MPR lines. 

Figure 3 shows the frequency shifts of the LO and 
TO inodes for different amounts of disorder under a 
compressive deformation e = —10%. The anticrossing 
is found to be smeared out as the disorder increases as 
found in the undeformed case^^. 

Fig. 3 shows also, that at a given magnetic field, the 
LO and TO frequency shifts are different as found in 
zero magnetic field^^. This feature is ascribed to the 
dependence of the electron-hole pair production on the 
strain, leading to long-lived (damped) phonons along 
the tensile (compressive) deformation axis^^. 

Several studies have given evidence of the electron con¬ 
centration dependence of the This depen¬ 

dence is shown hgure 4 where we plot the frequency shifts 
and the broadening of the LO and TO modes for different 

Fermi energy value. We denote ~ 1.22 hujo 

the Fermi energy at which the n = 0 LL is totally filled 
in unstrained graphene and (l -|- ^/2^ huQ ~ 

SAhuo corresponds to the totally occupied n = 1 LL®. 

As shown by Fig. 4, one can distinguish three regimes 
for the MPR behavior, depending on the Fermi energy: 

(i) for 0 < Ep < where both 0 <=l ±1 

transition are active, the anticrossing is well defined with 
a two strain-induced split lines. These lines are reminis¬ 
cent of the Raman spectrum obtained by Kashuba and 



FIG. 3. Frequency shifts of the LO and TO modes as a func¬ 
tion of the magnetic energy Aub in the undoped graphene at 
e = —10% and for different disorder amounts. The shifts and 
the broadenings of the LO modes, which is along the strain 
direction, are more pronounced than the corresponding TO 
modes due to lattice softening along the deformation axis. 

Fal’ko^^ at low carrier concentration and for linearly po¬ 
larized light. 

(ii) For < Ep < eP ~ SAhujo, the — 1 —!► 0 tran¬ 
sition is blocked by Pauli principle and only the 0 —>■ 1 
transition is active. 

The phonon broadening, in this doping regime, is 
strongly reduced compared to the case of low doping 
when the n = 0 is not totally filled. This reflects the 
reduction of the number of electron-hole pairs interact¬ 
ing with phonons. 

(iii) For Ep > eP, the n = I LL is totally occupied 
and the inter-LL transition 0 —> 1 is now stopped. The 
corresponding MPR disappears. We obtain a single line 
at the bare phonon frequency ujq corresponding to the 
undeformed lattice. The strain-induced splitting of the 
phonon line disappears since the latter is due to the 
coupling of phonons to electron-hole pairs. Such cou¬ 
pling can no more take place since the pair production 
is blocked by Pauli principle. Our results are consistent 
with those obtained by Kashuba and Fal’ko^^ in the 
high doping limit where two field independent Raman 
lines were reported for circular polarized light and a 
single line for linearly polarized light. The latter shows 
a relative shift compared to the bare phonon frequency 
wq. This shift, which does not appear in our results, is 
due to the shear strain component considered in Ref. 21. 
This component gives rise to a degeneracy lifting of the 
LO and TO phonon modes even in the absence of the 
coupling of phonon to electron-hole pairs^®. 

It is worth mentioning that the interaction of op- 
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FIG. 4. Frequency shifts and broadenings of the TO and LO 
modes as a function of the magnetic energy Hlub at e = 10% 
for different doping levels and for a disorder amount of 5 = 
0.02/kuo. The LO mode is along the strain axis. 


tical phonon with non-relativistic electrons in 2D sys¬ 
tems under a transverse magnetic field has attracted 
considerable interesF^^. This interaction induces an 
avoided crossing of the LL when the cyclotron reso¬ 
nance frequency matches an appropriate optical phonon 
frequency^®. This effect is substantially different from 
that studied in the present work where we consider the 
role of electron-phonon coupling on the phonon proper¬ 
ties and not on the electronic spectrum. 


B. MPR in uniaxial strained graphene under an 
inplane electric field 

To avoid the collapse of LL, the effective tilt param¬ 
eter W{(£’) should satisfy w^{E) < l24.25^ Given the 
expression of 'w^{E) (Eq. 13), the latter condition can be 
expressed as 

~ |0.6e-^-^| < 1 (32) 

where we considered the approximate value of ~ 0 . 6 e 
and put Wx hvp- 

Regarding the small values of the strain amplitude, 
the condition on the effective tilt parameter w^{E) < 1 
reduces to —^ < 1 as found by Lukose et al?'^. Since 
the resonance corresponding to 0 ^ ±1 transitions take 
place at i? ^ 30 T and given the graphene average Fermi 
velocity vp ~ 10 ®m.s“^, the electric field amplitude 
should not exceed 10 ^ V-m“^. 

Figure 5 shows the frequency shifts and the broad¬ 
enings of the LO mode as a function of the magnetic 
field energy for different strain values in the undoped 
graphene and under an electric field of 5 x 10® V.m“^. 
The shift line exhibits a double peak structure which 
gets more pronounced as the strain amplitude increases. 
According to our numerical results, no double peak 
is expected at a tensile (compressive) deformation 
below e = 15% (|e| = 10%) for a disorder amount of 
5 = 0.02 hujo- 

Figure 5 shows that, for the chosen disorder amount 
and doping, the MPR lines exhibit a remarkable behavior 
above the critical strain value. A double peak structure 
emerges as a fingerprint of the twofold valley degeneracy 
lifting induced by the effective tilt parameter w^{E). As 
discussed above, this valley degeneracy lifting is due to a 
different inter-LL spacing in both valleys induced by the 
inplane electric field (Eq. 15). Therefore, The inter-LL 
transitions 0 ^ ±1 take place at different resonant 
magnetic fields giving rise to the double resonant MPR 
lines. 

It is worth noting that, in the undeformed case 
(solid line in Fig. 5 ), there is no electric field induced 
double peak. The effective tilt parameter reduces, in 
the absence of strain, to w{E) = (Eq. 32) which is 
independent of the valley index. Therefore, the tilt of 
Dirac cones induced by the uniaxial strain in graphene 
is a substantial ingredient to lift the twofold valley 
degeneracy in the presence of the inplane electric field. 
The double peak structure of the MPR line emerges as 
a consequence of the simultaneous action of the Dirac 
cone tilt and the electric field. 

Figure 5 shows also that, the double peak effect is 
more pronounced for a compressive deformation. This 
can be understood from the resonance condition, given 
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FIG. 5. Frequency shifts (a) and broadenings (b) of the LO 
mode as a function of the magnetic energy fens for differ¬ 
ent strain strength in the undoped case for a disorder amount 
= 0.02 and under a uniform electric field of 5 x 10® V.m“^. 

TUJQ 

The LO mode is along the strain axis. The double peak struc¬ 
ture is due to the lifting of the twofold valley degeneracy. 


by Eq. 31, which can be written, in the presence of the 
inplane electric field, as: 

fujo = nhj% ^1 - (l - wl {E)) ^ (33) 

where we considered the approximate expressions of 
and Wy (Eq. 5). The effective tilt parameter Eq. 13 can 
be written as 


( 34 ) 

According to Eq. 34, the larger effective tilt parameter 
is obtained for a compressive deformation (e < 0) and in 
the ^ = — valley. The corresponding resonance is then 
expected to appear before that ascribed to the ^ = + 
valley. The difference between the valley resonance fields 
for compressive strain is larger than that obtained for a 
tensile deformation. The double peak structure is then 
expected to be more pronounced for compressive strain 
as shown in Fig. 5. 

As discussed in the previous section, the fine structure 
of the MPR depends on the disorder amount. The 
latter is then expected to be a key parameter for the 
occurrence of the double peak in the MPR line. Figure 
6 shows the frequency shifts and the broadenings of the 
LO phonon mode around the anticrossing for different 
disorder amount and for an electric field of 5 x 10® V-m“^. 

As shown by Fig. 6 (b), the double peak structure is 
smeared out by increasing the disorder. This effect could 
be understood from the expression of the phonon self¬ 
energy Eq. 19. The resonance condition reflecting the 
divergence of the self-energy depends on the disorder pa¬ 
rameter 6. The latter may, if increased, erase the anti¬ 
crossing and the relative double peak structure. 

The double peak separation Au}{B) is given by = 

—i-, where a;|j (Eq. 29) yields to the singular 
behavior of the self-energy given by Eq. 19. 

From this equation, one can deduce a simple condition for 
the occurrence of the double peak structure as a function 
of the disorder. Such condition could be written as: 

Auj{B) S 

Wo Wo 

The double peak is not smeared out by disorder, as far 
as this condition is satisfied. 

The keystone parameter for the occurrence of the 
double peak structure is the effective tilt of Dirac cones 
induced by the electric field. The latter should obey 
to the condition < 1 to avoid the collapse of the 
LL^^. Figure 7 shows the MPR double peak structure 
for different electric field amplitude satisfying the 
above mentioned condition. The double peak structure 
becomes more pronounced as the electric field increases. 
This behavior results from the electric field dependence 
of the effective tilt parameter w^(E) which is enhanced 
by increasing the electric field. The higher the electric 
field, the larger the difference between LL spacings in 
the valleys. 
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(b) 

FIG. 6. Frequency shifts (a) and the broadening's (b) of the 
LO phonon mode as a function of the magnetic energy Hlub 
for compressive deformation of —18% under a uniform electric 
field of 5 X 10® V-m“^ and for different disorder amounts. 


IV. CONCLUDING REMARKS 

We studied the F point longitudinal (LO) and trans¬ 
verse (TO) optical phonon modes, coupled to the inter- 
LL transition 0 ±1, around the magneto-phonon res¬ 

onance field i? ^ 30 T in uniaxial strained graphene. 
We derived the frequency shifts and the broadenings of 
the magneto-phonon modes taking into account the cou¬ 
pling of phonons to the strain modified electronic spec¬ 
trum. The latter is characterized by tilted Dirac cones 
■with anisotropic velocities. We found that, disregard¬ 
ing the shear strain component and the strain effect on 




Magnetic energy (Units of RcOq) 

(b) 

FIG. 7. Frequency shifts (a) and the broadenings (b) of the 
LO phonon mode as a function of the magnetic energy hujs 
for a compressive deformation of —20% under different electric 
field amplitudes and for disorder amount of S/huuo = 0.02. 


the phonon dispersion, the phonon-electron coupling con¬ 
tributes to the split of the MPR lines and to the shift of 
the resonant magnetic field. We also show that the criti¬ 
cal doping, and the corresponding filling factor, responsi¬ 
ble of the quench of an inter-LL transition are strain de¬ 
pendent. This effect may be used to measure the strain 
amplitude in the sample. Moreover, we found that, un¬ 
der an inplane electric field, the simultaneous effect of 
the Dirac cone tilt and the electric field, gives rise to a 
double peak structure in the MPR spectrum of an opti¬ 
cal phonon mode ^ (p, = LO, TO). This new structure 
is the signature of the valley degeneracy lifting resulting 
from different couplings of phonon to electron-hole pairs 
in the two valleys, where Landau level spacing are differ¬ 
ent. The MPR spectrum of graphene under a uniaxial 
strain and in the presence of an inplane uniform electric 
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field, is then expected to show four peaks, two for each 
phonon mode. We show that the double peak structure 
is substantially dependent on the disorder amount in the 
lattice, on the doping and on the electric field amplitude. 
We propose that this structure of the MPR spectrum cor¬ 
responding to the 0 ^ ±1 transitions, could be observed 
in pristine graphene for a uniaxial compressive deforma¬ 
tion of ~ 15% and under an inplane electric field of 5 x 10® 
V-m“^ around the resonant field B ^ 30 T. 
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